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We have observed the superfluid phase transition in a strongly interacting Fermi gas via high- 
precision measurements of the local compressibility, density and pressure down to near-zero entropy. 
Our data completely determine the universal thermodynamics of strongly interacting fermions with- 
out any fit or external thermometer. The onset of superfluidity is observed in the compressibility, 
the chemical potential, the entropy, and the heat capacity. In particular, the heat capacity displays 
a characteristic lambda-like feature at the critical temperature of T c /Tf = 0.167(13). This is the 
first clear thermodynamic signature of the superfluid transition in a spin-balanced atomic Fermi gas. 
Our measurements provide a benchmark for many-body theories on strongly interacting fermions, 
relevant for problems ranging from high-temperature superconductivity to the equation of state of 
neutron stars. 



Phase transitions are ubiquitous in Nature: water 
freezes into ice, electron spins suddenly align as materi- 
als turn into magnets, and metals become superconduct- 
ing. A pervasive feature of continuous phase transitions 
is their critical behavior, namely singularities in thermo- 
dynamic quantities: the magnetic susceptibility diverges 
at a ferromagnetic transition, the specific heat shows a 
jump at superconducting transitions pQ as well as at the 
superfluid transition of 3 He j2]. In 4 He, at the famous A- 
transition into the the superfluid state, the jump is even 
resolved, in zero gravity, to be a near-diverging, singular 
peak [3]. A novel form of superfluidity has been real- 
ized in trapped, ultracold atomic gases of strongly in- 
teracting fermions, particles with half- integer spin [IHB]. 
Thanks to an exquisite control over relevant system pa- 
rameters, these gases have recently emerged as a versatile 
system well suited to solve open problems in many-body 
physics [§]. However, while superfluidity has been estab- 
lished via the observation of vortex lattices in rotating 
gases [7J, no clear thermodynamic signature of the su- 
perfluid transition has previously been observed. 

Initial measurements on the thermodynamics of 
strongly interacting Fermi gases have focused on trap av- 
eraged quantities [8HTU] in which the superfluid transition 
is inherently difficult to observe. It is also challenging to 
reveal the critical behavior through the study of local 
thermodynamic quantities. The emergence of the con- 
densate of fermion pairs in a spin-balanced Fermi gas is 
accompanied by only minute changes in the density 0]. 
Therefore, quantities that involve integration of the den- 
sity over the local potential, such as the energy E [TT] 
and the pressure P |12j . are only weakly sensitive to the 
sudden variations in the thermodynamics of the gas that 
one expects near the superfluid phase transition [TSJ . 

A thermodynamic quantity involving the second 
derivative of the pressure P is expected to become sin- 
gular at the second order phase transition into the su- 
perfluid state. An example is the (isothermal) compress- 



ibility k = ^fp|r, the relative change of the gas density 
n due to a change in the pressure P. As the change in 
pressure is related to the change in chemical potential 
/j, of the gas via dP = ndfi at constant temperature, 
K = ^2 §tJt|t is a second derivative of the pressure, and 
thus should reveal a clear signature of the superfluid tran- 
sition. 

In the present work, we report on a high-precision mea- 
surement of the local compressibility, density and pres- 
sure across the superfluid phase transition in the unitary 
Fermi gas, realized using a trapped, two-spin mixture 
of fermionic 6 Li atoms at a Feshbach resonance [H [6]. 
The combination of these three directly measurable, lo- 
cal quantities determines the entire thermodynamics of 
the homogeneous gas as we show below. Thermometry 
of strongly interacting Fermi gases was known to be no- 
toriously difficult. Our method solves this problem. It 
is general and applies to other systems as well, such as 
Bose gases in three or two dimensions. 

In order to determine the thermodynamic properties 
of a given substance, the general goal is to measure an 
equation of state (EoS), such as the pressure P(/jI, T) as a 
function of the chemical potential fi and the temperature 
T. Equivalently, replacing the pressure by the density 
n = §^|tj one can determine the density EoS n(fj,,T). 

In our experiment we directly measure the local gas 
density n(V) as a function of the local potential V from 
in-situ absorption images of a trapped 6 Li gas. The trap- 
ping potential is cylindrically symmetric, with harmonic 
confinement along the axial direction. This symmetry 
allows to find the 3D density via the inverse Abel trans- 
form of the measured column density [14] . Other than 
cylindrical symmetry, no other assumption on the shape 
of the potential is made. Instead, the local potential is 
directly measured via the atomic density distribution and 
the accurately known harmonic potential along the axial 
direction of the atom trap. 

The compressibility k follows as the change of the den- 
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sity n with respect to the local potential V experienced 
by the trapped gas. The change in the local chemical 
potential is given by the negative change in the local 
potential, d/tt = — dV, and hence the compressibility is 

„ _ 1 An I 

K - -^av\ T - 

The compressibility allows replacing the unknown 
chemical potential /j, in the density EoS n(/z, T) by its 
known variation in the atom trap, yielding the equation 
of state n(K, T). In lieu of the a priori unknown tempera- 
ture, we determine the pressure P(V) — d//n(//) = 
Jy dVn(V) that is simply the integral of the density 
profile over the potential [121 HE]- The resulting equa- 
tion of state n{n, P) contains only quantities that can 
be directly obtained from the density distribution. This 
represents a crucial advance over previous methods that 
require the input of additional thermodynamic quanti- 
ties, such as the temperature T and the chemical poten- 
tial fi, whose determination requires the use of a fitting 
procedure or an external thermometer, as in [111 112] . 

We normalize the compressibility and the pressure by 
the respective quantities at the same local density for a 
non-interacting Fermi gas at T = 0, kq = | ^gr- and 

Pq = ynEp, where Ep = - ^ 3 Z — is the Fermi en- 
ergy and m is the particle mass, yielding k = k/k and 
p = P/Pq- In general, k could depend on dimensionless 
quantities other than P/Pq, such as the interaction pa- 
rameter na 3 of the gas, where a is the scattering length. 
However, for dilute gases at the Feshbach resonance, the 
scattering length diverges and is no longer a relevant 
length scale. In the absence of an interaction-dependent 
length scale, the thermodynamics of such resonant gases 
are universal |17j , and k is a universal function of p only. 
Every experimental profile n(V), whatever the trapping 
potential, the total number of atoms or the temperature, 
must produce the same universal curve k versus p. By 
averaging many profiles, one obtains a low-noise deter- 
mination of k(p). 

Our method has been tested on the non-interacting 
Fermi gas that can be studied in two independent ways: 
In spin-balanced gases near the zero-crossing of the 
scattering length and in the wings of highly imbal- 
anced clouds at unitarity, where only one spin state is 
present locally. Both determinations yield the same non- 
interacting compressibility EoS. 

Fig. 1 shows the resulting compressibility equation of 
state k(p). In the high-temperature and, equivalently, 
low-pressure regime, the pressure (and hence all other 
thermodynamic quantities) allows for a Virial expansion 
in terms of the fugacity [T3]: P/3A 3 = 2V\ bje j ^, 
with the nth order Virial coefficients b n . It is known 
that b x = 1, b 2 = 3V2/8 [20], and 63 = -0.29095295 \19\. 
Our data agrees excellently with the Virial expansion. 
Fixing b 2 and 63, our measurement yields a prediction for 
64 = +0.065(10), in agreement with [T2], but contradict- 
ing a recent four-body calculation that gives a negative 
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FIG. 1. Compressibility k versus pressure P of the uni- 
tary Fermi gas. Red solid dots show k and P normalized by 
their respective values ko and Po for a non-interacting Fermi 
gas at the same density. The error bars show one standard 
deviation. The blue solid line shows the 3rd order Virial ex- 
pansion. Black open squares (black solid line) denote data 
(theory) for a non-interacting Fermi gas. The black dashed 
line shows the relation ft = 1/p that must be obeyed at zero- 
temperature both for the non-interacting gas (ft = 1/p = 1) 
and the unitary gas (ft = 1/p = l/£) (dotted lines). A gray 
band marks the uncertainty region for the T — value of 
ft = l/£ and p = £. The blue dashed curve is a model for the 
EoS of the unitary Fermi gas (above T c : interpolation from 
the Monte-Carlo calculation |18| . below T c : BCS theory in- 
cluding phonon and pair- breaking excitations). Green solid 
line: Effect of 2/j.m optical resolution on the model EoS. 

sign IS]. 

At degenerate temperatures, the normalized compress- 
ibility rises beyond that of a non-interacting Fermi gas, 
as expected for an attractively interacting gas. A sudden 
rise of the compressibility at around p = 0.55, followed 
by a decrease at lower temperatures marks the superfluid 
transition. Indeed, the theory of superfluid phase tran- 
sitions in three dimensions implies a singularity of the 
compressibility at the transition. This would appear as a 
sudden change in slope of n(V) that is rounded off by the 
finite resolution of our imaging system. Below the tran- 
sition temperature, the decrease of the compressibility 
is consistent with the expectation from Bardccn-Cooper- 
Schrieffer (BCS) theory where single-particle excitations 
freeze out and pairs form (see model in Fig. 1). 

As T — > 0, the Fermi energy Ep is the only intensive 
energy scale, and so the chemical potential must be re- 
lated to Ep by a universal number, [i = £,Ep, where £ 
is known as the Bertsch parameter [S] E] . It follows that 
at T = 0, k = 1/p = l/£. The extrapolation of the 
low-temperature experimental data for k(p) towards the 
curve k — 1/p gives £ ~ 0.37, a value that we find consis- 
tently for the normalized chemical potential, energy and 
free energy at our lowest temperatures (see below). 

From the universal function k{p) we obtain all other 
thermodynamic quantities of the unitary gas. First of all, 
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FIG. 2. A Normalized compressibility k — ^kuEf and B specific heat per particle Cv/Nks of a unitary Fermi gas versus 
reduced temperature T/Tf, shown in solid red circles. The black solid curve shows the theory for a non-interacting Fermi 
gas, the blue solid curve shows the third order Virial expansion for the unitary gas. Black open squares denote data for a 
non-interacting Fermi gas. The green solid line shows the model from Fig. 1, again including 2/im imaging resolution. C 
Condensate fraction at unitarity as determined from a rapid ramp to the molecular side of the Feshbach resonance. The onset 
of condensation coincides with the sudden rise of the specific heat. All error bars show one standard deviation. 



it allows translation between the pressure thermometer 
p and the more familiar normalized temperature T/Tp 
(where k B Tp = Ep). Considering that the change in 
p with T/Tp at constant temperature is related to the 
change in pressure with density and thus the compress- 
ibility, one finds d (7y^ F ) = f^ 1 {p — A) and thus by in- 



tegration [T3] 
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Here, (T/Tp)i is the normalized temperature at an initial 
normalized pressure pi that can be chosen to lie in the 
Virial regime validated above. 

Thanks to the relation E = f-PV valid at unitarity [T7] . 
we can also directly obtain the heat capacity per particle 
at constant volume V 1131. 
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Fig. 2 shows the normalized compressibility and the 
specific heat as a function of T/Tp. At high tem- 



peratures, the specific heat approaches that of a non- 
interacting Fermi gas and eventually Cy = ^Nk B , the 
value for a Boltzmann gas. A dramatic rise is observed 
around T C /T F 0.16, followed by a steep drop at lower 
temperatures. Such a A-shaped feature in the specific 
heat is characteristic of second order phase transitions, 
as in the famous A-transition in 4 He |23j . Jumps in 
the specific heat are also well-known from superconduc- 
tors pQ and 3 He [2], To our knowledge, this is the 
first time that a specific heat jump has been directly 
measured in an ultracold atomic gas. Previously, such 
jumps had only been inferred from derivatives to fit func- 
tions that implied a jump [TUl El]- We do not expect 
to resolve the critical behavior very close to T c , given 
by C v = C ns + A +/ ^\T - T c \~ a with C ns the non- 
singular part of the specific heat, the critical exponent 
a —0.012 and amplitudes A + /_ above and below the 
transition [25]. In our trapped sample the critical region 
is confined to a narrow shell due to the spatially vary- 
ing chemical potential. Based on the estimate in |26j . 
the thickness of the critical shell is 1% of the cloud size. 
The finite resolution of our imaging system (2 /im or 
about 5% of the cloud size in the radial direction) suf- 
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FIG. 3. A Chemical potential fi, energy E and free energy F of the unitary Fermi gas versus T/Tp. i-i is normalized by 
the Fermi energy Ep (red solid circle), E (black solid circle) and F (green solid circle) are normalized by Eo = ^NEp. At 
high temperatures, all quantities approximately track those for a non-interacting Fermi gas, shifted by — 1 (dashed lines). 
The peak in the chemical potential signals the onset of superfluidity. In the deeply superfluid regime at low temperatures, 
[i/Ep, E/Eo and F/Fo all approach £. B Entropy per particle. At high temperatures, the entropy closely tracks that of a 
non-interacting Fermi gas (solid line). Below T/Tp = 0.1, the entropy per particle reaches values < 0.04 fcg. The open squares 
are from the self-consistent T-matrix calculation \22 . A few representative error bars are shown, representing one standard 
deviation. 



fices to explain the rounding of the singularity expected 
from criticality. The rounding reduces the jump in the 
heat capacity at the transition. From our data, we ob- 
tain AC/N = C s /N - C n /N, where C s /N {C n /N) is 
the specific heat per particle at the peak (the onset of 
the sudden rise). Due to the finte resolution, our mea- 
surement sets a lower bound on the jump in the specific 
heat, AC/C n > l.Oj^. Considering the strong interac- 
tions, this is surprisingly close to the BCS value of 1.43. 
Below T c , the specific heat is expected to decrease as 
~ exp(— Ao/fceT) due to the pairing gap A . At low 
temperatures T T c the phonon contribution oc T 3 
dominates [22J. This behavior is consistent with our data, 
but the phonon regime is not resolved. 

To validate our in-situ measurements of the super- 
fluid phase transition, we have employed the rapid ramp 
method to detect fermion pair condensation [37J HH]- 
The results are shown in Fig. 2C, and demonstrate 
that the onset of condensation and the sudden rise in 
specific heat and compressibility all occur at the same 
critical temperature, within the error bars. Previous 
experimental determinations of T c /Tp of the homoge- 
neous unitary Fermi gas relied on a temperature cali- 
bration that did not agree with the Virial expansion |llj , 
or obtained Tp from the /i-derivative of a pressure fit, 
and T c from a construction implying a first-order tran- 
sition [12] . Here we determine T c /Tp directly from the 
density profiles, finding a sudden rise in the specific heat 



and the onset of condensation at T c /Tp = 0.167(13). 
This value is determined as the midpoint of the sudden 
rise, and the error is assessed as the shift due to the un- 
certainty of the Feshbach resonance [13] • This is in very 
good agreement with theoretical determinations, for ex- 
ample, the self-consistent T-matrix approach that gives 
T c /Tp w 0.16 [221 [29], and Monte-Carlo calculations 
which give T C /T F = 0.173(6) [30] and 0.152(7) [31]. It 
disagrees with T C /T F = 0.23(2) [32], but is close to a later 
determination by the same group of T c /Tp < 0.15(1) [53] , 
There is a current debate on the possibility of a pseudo- 
gap phase of preformed pairs above T c [T2J [34] . A pair- 
ing gap for single-particle excitations above the transition 
should be signaled by a suppression of the specific heat 
already above T c , which is not observed in our measure- 
ments. 

From the definition of the compressibility k = ^2 §^|t 
we can obtain the reduced chemical potential fi/Ep as 
a function of the reduced temperature, see Figure 3A 
and [T3]. Here, for the first time, this function is obtained 
from measured quantities, rather than from numerical 
derivatives of data that involved uncontrolled thermom- 
etry pi]- Around T /Tp ~ 0.25 to 1, the chemical po- 
tential is close to that of a non-interacting Fermi gas, 
shifted by (£„ — l)Ep due to interactions present already 
in the normal state, with £„ « 0.45. However, unlike a 
normal Fermi gas, the chemical potential attains a max- 
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FIG. 4. A Density and B pressure of a unitary Fermi gas versus fi/ksT, normalized by the density and pressure of a non- 
interacting Fermi gas at the same chemical potential [i and temperature T. Red solid circles: experimental EoS. Dashed lines: 
low-temperature behavior with £ = 0.364, 0.376 and 0.389. Black dashed line: low-temperature behavior from the £ upper 
bound £ = 0.383 [35]. Green open circles and black dashed line at 1.0: MIT experimental density and pressure, and theory for 
the ideal Fermi gas. Blue solid squares (blue band): Diagrammatic Monte Carlo 18 for density (pressure). Solid green line: 
3rd order Virial expansion. Open black squares: self-consistent T-matrix |22| . Open green circles: [32j. Orange star: [30]. Blue 
star: Solid diamonds: ENS experiment Q2]. Open diamonds: Tokyo experiment [11] , 



imum of fi/E F = 0.42(1) at T/T F = 0.171(10), and then 
decreases at lower temperatures. This is expected for a 
superfiuid of paired fermions |22j . As the temperature 
is increased from zero in a superfiuid, phonons (sound 
excitations) emerge that increase the chemical potential 
(i. In addition, fermion pairs start to break and single 
fermions contribute increasingly to the chemical poten- 
tial with increasing temperature. At T c , n/E F must have 
a sharp change in slope, as d(/j,/ E F )/d(T/T F ) involves 
the singular compressibility. Indeed, the self-consistent 
T-Matrix calculation shows a very clear peak in n/E F 
near T c |22j . in agreement with our observation. At low 
temperatures, the reduced chemical potential fi/E F sat- 
urates to the universal value £. As the internal energy 
E and the free energy F satisfy E(T) > E(T = 0) = 
lN£E F = F(T = 0) > F(T) for all T, the reduced 
quantities f E = = p and f F = = |^ - §p 

provide upper and lower bounds for £ |36j , shown in Fig. 
3A. Taking the coldest points of these three curves and 
including the systematic error due to the effective inter- 
action range, we find £ = 0.376(5). The uncertainty in 
the Feshbach resonance is expected to shift £ by at most 
2% [13]. This value is consistent with a recent upper 
bound £ < 0.383 [35], is close to £ = 0.36(1) from a 
self-consistent T-matrix calculation [22] , and agrees with 
£ = 0.367(9) from an epsilon expansion [37]. It lies be- 
low earlier estimates £ = 0.44(2) [38 and \ = 0.42(1) [39] 
via fixed-node quantum Monte-Carlo that provide upper 



bounds on £. Our measurement agrees with several less 
accurate experimental determinations 5 , but disagrees 
with the most recent experimental value 0.415(10) that 
was used to calibrate the pressure in |12j . shown in Fig. 
4B. 

From the energy, pressure and chemical potential, we 
can obtain the entropy S = ^(E + PV — fJ.N). Shown 
in Figure 3B is the entropy per particle S/Nks — 
"t~(p ~ Bp) as a function of T/T F . At high tempera- 
tures, S is close to the entropy of an ideal Fermi gas at 
the same T/T F . Down to T c , neither the non-interacting 
nor the unitary Fermi gas has S/N -C fcs- Also, the 
specific heat Cy is not linear in T. Thus it is question- 
able to identify the normal regime as a Landau Fermi 
Liquid, although some thermodynamic quantities agree 
surprisingly well with the expectation for a Fermi liquid 
(see [12] and [13]). Below about T/T F = 0.17 the en- 
tropy starts to strongly fall off compared to that of a 
non-interacting Fermi gas, which we again interpret as 
the freezing out of single-particle excitations due to for- 
mation of fermion pairs. Far below the critical tempera- 
ture for superfluidity, phonons dominate. They only have 
a minute contribution to the entropy [52], less than 0.02 
is at T/T F = 0.1, consistent with our measurements. 
At the critical point we obtain S c — 0.73(13)Afcs, in 
agreement with [22]. It is encouraging for cold atom ex- 
periments that we obtain very low entropies, less than 
0.04iVfcs, far below critical entropies required to reach 
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magnetically ordered phases of fermions in optical lat- 
tices. 

From the chemical potential fi/Ep and T/Tp = 
' we finally obtain the density EoS 
with the de Broglie wavelength 
The pressure EoS follows as P(/x, T) = 



(3tt 2 ) 2 / 3 (nA 3 ) 2 / 3 ' 

n(/i,T) = Jy/„(/3/i) 

2ttR 2 



A = 



%k B T ' 



^/ P (/3/i) with / P - l^fpfnifo). Fig. 4 shows the 
density and pressure normalized by their non-interacting 
counterparts at the same chemical potential and tem- 
perature. For the normal state, a concurrent theoret- 
ical calculation employing a new Monte-Carlo method 
agrees excellently with our data JT5j. Our data deviates 
from a previous experimental determination of the pres- 
sure EoS |12j that was calibrated with an independently 
measured value of £ = 0.415(10) [40] - a value which 
agrees neither with our value nor a recent theoretical up- 
per bound 35] . Our data also disagrees with the energy 
measurement in |llj that used a thermometry inconsis- 
tent with the Virial expansion |10) . Around the critical 
point the density shows a strong variation, while the pres- 
sure, the integral of the density over \x at constant T, is 
naturally less sensitive to the superfluid transition. 

In conclusion, we have performed a high-precision ex- 
perimental study of the thermodynamics of a unitary 
Fermi gas across the superfluid transition. For the first 
time, we observe direct signatures of the superfluid phase 
transition in various thermodynamic quantities, such as 
a sudden rise in the compressibility and the specific heat. 
The critical temperature for superfluidity at unitarity is 
found to be T c = 0.167(13)Tp, the critical entropy is 
0.73(13)-/Vfcs, and we set a lower bound on the jump in 
the specific heat, AC/C n > 1.0^. We also find a new ex- 
perimental value for the Bertsch parameter £ = 0.376(5), 
subject to at most a 2% shift due to the uncertainty in 
the Feshbach resonance position. We have demonstrated 
that precision many-body determinations of phase tran- 
sitions involving fcrmionic atoms are possible at the few 
percent level and without any fits or input from theory, 
enabling validation of theories for strongly interacting 
matter. Similar unbiased methods can be applied to 
other systems, for example, two-dimensional Bose and 
Fermi gases or fermions in optical lattices. 
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SUPPLEMENTAL MATERIALS AND METHODS 

Obtaining density versus potential curves from 
absorption images of trapped samples 

In this section we provide details on the determination 
of the density from in-situ images of trapped samples. 
The experimental set up is described in [Slj . Fermionic 
6 Li is cooled to degeneracy via sympathetic cooling with 
23 Na. A two-state mixture of the two lowest hyperfine 
states of 6 Li is prepared and brought into the strongly- 
interacting regime at a broad Feshbach resonance at 
834.15 G in a hybrid magnetic and optical trap, where 
the mixture is further cooled. We measure the column 
density of the gas via high-resolution in situ absorption 
imaging. 

The trapping potential is provided by a magnetic field 
saddle potential that is confining along the z-direction, 
and a gaussian laser beam (wavelength 1064 nm, waist 
w = 120 (im) propagating along z. Due to the large 
waist and correspondingly large Rayleigh range of the 
gaussian beam, the trapping potential is to an excellent 
approximation harmonic in the z-direction, with a mea- 
sured trapping frequency of uj z — 2n ■ 22.83(5) Hz. Waist 
w and power V of the laser beam determine the trapping 
potential along the radial direction. The total potential 
can be excellently modeled as 

V(p,z) = V p (p) + V z (z) 

_ 1 2 2 1 2 2 r -2p 2 /w 2 lo\ 

= -muj z z — -muj,p —a- — H/ , (3) 
2 z 4 z \w 2 v ; 

where a is the polarizability. We can safely neglect in 
this model the 0.01% contribution to the axial trapping 
frequency from the variation of the laser power along 
the z-direction. We do not rely on the waist and power 
measurement of the laser beam but in fact measure the 
trapping potential directly using the atomic distribution 
itself, as will be outlined below. 

To obtain the 3D density from the measured optical 
density OD(x,z) = oniT)(x,z) with effective absorption 
cross section cr, we employ the inverse Abel transform for 



each z [32]: 

n 3 D{x,z) = -- dp — — (p,z) 

n Jx VP _ X " P 

Cylindrical symmetry, i.e. the equality of the waist of 
the optical potential along the transverse x and y di- 
rection, is measured via a beam profiler to be better 
than 1% and confirmed via the observation of long-lived 
vortices |S3j . Cylindrical symmetry is not required: As 
long as equipotential lines in the (x,y) plane transverse 
to the z-direction can be parameterized as an ellipse 
x 2 /a 2 + y 2 /b 2 = const, the inverse Abel formula for the 
3D density would be modified only through an overall 
factor b/a: 

b 7T J x y/ p 2 _ x 2 dp 

The absolute scale of density is calibrated through nor- 
malization via the measured effective absorption cross 
section of a non-interacting Fermi gas, so this prefac- 
tor, even if different from unity, is simply absorbed in 
the effective absorption cross section. Absorption imag- 
ing is performed at low intensity ///sat < 0.07 to avoid 
saturation, and the total number of photons scattered 
is kept small (around 4) by using a short pulse dura- 
tion to avoid optical pumping and an increasing Doppler 
shift |S4] . Simulations have confirmed that non- linear ef- 
fects in our imaging procedure distort the density profile 
by much less than 1%. 

To obtain a direct measurement of the trapping po- 
tential, one can employ the local density approxima- 
tion: equidensity lines must be equipotential lines, and 
the potential is excellently known in the axial (z) di- 
rection, thus calibrating the potential everywhere. For 
this measurement, one may average many independent 
3D density profiles, even if the atom number or tem- 
perature vary from shot to shot. The averaged density 
along z, ra^f'(0, z) — n^'(V(0, z)) can be inverted to 
find V(ti3d), so that the density at any point (p, z) in 
the trap gives V(p, z). The resulting potential can be fit 
with our model, and results in a waist that agrees with 
the measured waist to within one percent, and a power 
that agrees with the measured power to within 10% (a 
typical error for standard power meters). Thus anhar- 
monicities in the trapping potential have been incorpo- 
rated into the experimental analysis of density profiles 
exactly. Previous measurements attempted to charac- 
terize anharmonic traps via measured trapping frequen- 
cies |S51 IS6j . However, anharmonicities cause such fre- 
quencies to depend on the trap filling, rendering that 
method less reliable. Knowing the trapping potential, a 
single experimental density profile n 3 D{p,z) can now be 
averaged over equipotential lines V(p,z) — Vo, yielding 
a low-noise determination of the density as a function of 
the local potential Vq. 



The local density approximation (LDA) is well fulfilled 
as long as the cloud size (typically ~ 40 fim in the narrow, 
radial direction) is large compared to the characteristic 
length of spatial correlations in the gas, on the order 
of ~ l/k F = 400 nm [ST. Violations of LDA due to 
criticality close to the superfluid transition have not been 
detected, likely due to our finite optical resolution (about 
2 — 3 /xm). 

The profiles n(V) for varying atom numbers, trap pa- 
rameters and degree of evaporative cooling are then used 
to calculate the normalized compressibility and the pres- 
sure, as described in the text. 



Compressibility equation of state and relation to the 
specific heat and reduced temperature 

The compressibility K, the pressure P and the density 
n are related by an equation of state K,(n,P). At uni- 
tarity, all these thermodynamic quantities are related to 
universal functions of /3/i. Let us define P = 

then n = ^f P (X) = ^/„(A), and nn 2 = ^ 
with X — ftfi. 

The experiment determines the density as a function 



of the local potential, n(V). We have 



dn 



On 
I) It 



— gy . We normalize k by the compressibility of a non- 
interacting, zero-temperature Fermi gas Kq = | . We 
dE^ = _9E fi where Ep _ H*(3™nf'° 



find k ee — 



k Ofl dV ' " iiVjiv - 2m 

the Fermi energy, with m being the mass of the atom. 

At unitarity, k is a universal function of T/Tp or, 
equivalently, of nX 3 . If we replace the temperature by 
the pressure P, another dimensionless quantity that we 
can form is p — , which is the pressure normalized 

by the zero-temperature limit for a non-interacting gas. 
One can trade the unknown thermometer T/Tp for the 
known, new thermometer p. For every profile, one can 



determine k(V) 



dE F (V) 



andp(V) 



/~ dV'n(V) 



and 



dV ™^F\y ) — Yn(V)E F (yy 
then plot k vs p. Even though obtaining k involves tak- 
ing a derivative of the data, accumulation of data reduces 
the noise. Data from all profiles within a bin of p is av- 
eraged to give the EoS. The statistical error bar is taken 
to be the standard error. 

and 



In terms of the function fp we have k = l^fr 

Jr 3 T f' p 



P 



| ~ 4r • The different thermometers, p. X 



Pfi and 



T/Tp depend on each other in the following way: 



dp 
dX 
d(T/T F ) 
dX 
dp 



5 T 
2Tp 



(1 - kp) 

T \ 2 
Tf) 



5T F 



1 



P 



(4) 



(5) 



(6) 



d{T/T F ) 

Note that apart from a prefactor, the last equation is just 
the heat capacity per particle (using E = |PV valid at 



unitarity) : 
C v 1 



dE 



k B N k B N dT 



dp 



N,V 



5 d (T/Tp) 



3 Tp 
2~T~ 



dp 



Further simplifying, we have j in(T /Tp ) 
which can be integrated to give 



Up 



T 
Tp~ 



dp 

P- 



(7) 
i \ 



(8) 



where k is known as a function of p, and (T/Tp)i is the 
temperature at the initial normalized pressure pi. This 
function relates the thermometer T/Tp to the pressure 
thermometer p. 

There are several ways to obtain vs T/Tp. One can 
consider function of T/Tp and obtain 



T/Tp 



Ti/Tp 



^ F k \ T 



(9) 



or one can consider T/Tp as a function of p and obtain 



/3/i = (/3/x)i 



Together with T/Tp 



dp 



1 



T 1 



Kp 



(10) 



(37T 2 ) 2 /3 f>2/ 



is f'ptfin) = WfJ,). We then get directly f P = 0^/', 



: , this now gives 

2 ~ Tp 



and fp = ^Tf^f'pk. In addition, for the entropy per 

particle we have S/Nks = ^f-(p— ^). 

We validate the method via the measurement of the 
equation of state of a non-interacting Fermi gas. This 
provides a check for the effective absorption cross section: 
If the cross section had been determined too large by a 
factor of 7, k oc j~ 2 / 3 would be too small, and p oc 7 2 / 3 
too large. Note, however, that the determination of T/Tp 
from k(p) is insensitive to errors in the absorption cross 
section, as all factors of 7 cancel in Eq. [8] Systematic 
errors that are unique to the unitary gas are discussed in 
the following. 



Systematic errors from non-universal scattering 
properties 

In the following we discuss systematic errors in our 
experiment from possible non-universal behavior. The s- 
wave scattering amplitude f(k) for atom-atom collisions 
can be written as 



1 



i(k) — ik 



where u(k) — —- + \r e k + a is the scattering 
length, and r e the effective range. The latter can gen- 
erally differ from the range of the interatomic potential, 
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FIG. SI. The effect of uncertainty in the Feshbach resonance 
position on the equation of state k versus p. The experimental 
result at 834.15 G is shown in red solid circles. This would be 
the EoS of the unitary Fermi gas if the Feshbach resonance is 
located at Bo — 834.15 G. If the true location of the resonance 
is at ±1.5 G away, one can estimate the true EoS for the 
unitary Fermi gas using the contact (see text). The result is 
shown as triangles pointing up (down) if the true resonance 
position is at 832.65 G (835.65 G). Also shown are the third 
order virial expansion (blue solid curve), the non-interacting 
EoS (black solid curve), and the curve k = 1/p on which the 
curve k versus p must end. 



the normalized chemical potential, energy and free en- 
ergy, pressure and entropy, which result from different 
combinations of density and pressure. All errors grow 
as the temperature is lowered, as the contact increases. 
At the lowest temperatures, the error in the density is 
3%, in the pressure 1%, in the chemical potential (and 
therefore £) 2%, in the normalized energy (free energy) 
4% (1%). As it turns out, the entropy per particle S/N 
is only very weakly sensitive to the uncertainty in the 
resonance position, as the systematic error in the pres- 
sure and density cancel to a high degree. At the lowest 
temperatures, the error has grown to only ±0.08fcs. At 
higher temperatures, the systematic error in S/N from 
the error in Bq is much smaller than the statistical error 
bars. Note that AB is much larger than the variation of 
the magnetic field along the axial direction of the cloud 
of about 15 mG due to the magnetic confinement, which 
contributes only a 0.01% error in £. 



Effective range correction 

The effective range near a Feshbach resonance at scat- 
tering length a can be modeled as [58] : 



b. To be in the universal regime at unitarity where the 
universal relation nA 3 = f(j3/i) holds, we require that 
k F a ^> 1, k F r e <C 1, k F b -C 1 [58]. Otherwise, the 
additional length scales a, r e and b will feature in the 
more complex equation of state that would be written 
as nA 3 = f(f3(i 7 \/a,\/r e ,\/b). We discuss the possible 
influence of the various terms separately in the following. 



Effect of the uncertainty in the Feshbach resonance position 

The Feshbach resonance between the lowest two hy- 
perfine states in 6 Li lies at (834.15 ± 1.5) G IS9 . The 
systematic error from the uncertainty in the Feshbach 
resonance position of AB = 1.5 G can be estimated us- 
ing the contact [SlOj . As the scattering length varies, the 
change in the energy of the system is given by 



dE 



da 1 J s N v 47rm 



CV 



where C is the contact density. From this relation, the 
change in the pressure and, via n = §^|t, the change 
in density with a can be derived. Diagrammatic Monte- 
Carlo calculations yield the temperature dependence of 
the contact [Sllj . We can thus predict the true density 
profiles starting with our measured normalized densities. 
The result is shown as the triangles pointing up and down 
in Fig. SI. 

From these upper and lower bounds on the true equa- 
tion of state we can deduce upper and lower bounds for 



r e = -2^(l-^) 2 + ^-- 
a \/7r a 



where J?* = 0.0269nm and b = 2.1 nm for 6 Li [55]. At 
our typical densities, 1/fcj? ~ 400 nm and thus k F R* < 
10~ 4 . As this parameter is very small, errors due to 
spatial variation of the scattering length, mediated by the 
optical trap, are vanishingly small |S8j . As a 3> b across 
the Feshbach resonance, the effective range is reduced to 
a constant r e — 4.7 nm for essentially all magnetic fields 
£> 650G [5T2] . 

At our highest densities at low temperatures, kpr e — 
0.012. All thermodynamic potentials and thus also the 
ground state energy E — ^£NEp are expected to de- 
pend linearly on this parameter (S 13] . From the Quan- 
tum Monte-Carlo values of £ versus effective range given 
in |S14] one can deduce a linear behavior of the upper 
bound on £: 

S(k F r e ) = C(r e = 0) + 0.26(6) • k F r e + ... 

With this dependence of £ on the effective range, the error 
on £ in our experiment is at most +0.003 or +0.8%. This 
is small compared to the error due to the uncertainty in 
the Feshbach resonance position. 



Energy vs entropy of the harmonically trapped gas 

In Fig. S2 we show the energy versus entropy for 
a harmonically trapped Fermi gas at unitarity, deter- 
mined from our experimental determination of the homo- 
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FIG. S2. Energy vs entropy of a harmonically trapped gas at 
unitarity, deduced from our measured bulk equation of state 
(red solid dots). The curve is smooth, as required by thermo- 
dynamics. Blue diamonds: Duke experiment jS15|, for which 
an interaction correction has been applied to the entropy. 

geneous density EoS and temperature. The data are com- 
pared with the earlier measurement of the same quanti- 
ties by the Duke group [SI 5], after applying a correction 
of the entropy for the finite interaction strength in that 
measurement. The agreement is excellent. Note that 
the energy vs entropy curve is smooth and on its own 
does not allow for a determination of the critical energy 
or entropy, in contrast to earlier interpretations of the 
Duke data |S15| . From the bulk thermometry, we obtain 
£ c ,t rap = 0.698(23)A^ jtrap and S Cjtrap = 1.70(10)JVfr B . 

Discussion on possible Fermi Liquid behavior 

In Fig. S3 and S4 we plot the density and pressure, 
normalized by that of a non-interacting Fermi gas at 
zero temperature, versus (T//ii) 2 . The linear behavior 
of the normalized density and pressure at high temper- 
atures is reminiscent of a Fermi liquid. However, in the 
range shown, T / fi is not a small parameter. The su- 
perfluid transition leads to a dramatic upturn at low T 
in the normalized density at low temperatures. This is 
directly related to the downturn in the chemical poten- 
tial seen in Fig. 3 A, as n(fj,,T)/no([M,0) oc (fj,/E F )~ 3 / 2 . 
The minimum value of the normalized density occurs at 
T/fi = 0.41(5), close to T c /n c = 0.40(3), where T c is 
determined from the midpoint in the specific heat jump 
and validated via the condensate fraction measurement. 

At low temperatures, the normalized density reaches 
the zero-temperature value l/£ 3 / 2 . The pressure does 
not show an upturn but smoothly attains a limiting value. 
The smooth behavior is expected for a second-order tran- 
sition, where first derivatives of the pressure with re- 
spect to the chemical potential are continuous across 
the transition. The intercept of a straight line at the 
limiting value P(fi,0)/Po(fi,0) = l/£ 3 ^ 2 with the linear 



4.8 




0.0 0.2 0.4 0.6 0.8 

(k B T/u) 2 

FIG. S3. Density n(fj,, T), normalized by the density no(/x, 0) 
of the non-interacting Fermi gas at zero temperature and same 
chemical potential fi, versus (T//i) 2 . The solid blue line de- 
notes the zero-temperature limit, n(/x, 0)/no(/z, 0) = l/f 3,/2 . 
The solid red line is a linear fit, resembling Fermi liquid be- 
havior of the density above T c . 
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FIG. S4. Pressure P(/z, T), normalized by the pressure 
Po(n,0) of the non-interacting Fermi gas at zero tempera- 
ture and same chemical potential fj,, versus (T//J,) 2 . Red solid 
circles are the experimental data from this work. Black solid 
circles are from [S6j . The solid blue line denotes the zero- 
temperature limit, P(/x, 0)/Po(a i ) 0) = l/£ 3 / 2 , the red solid 
line is a linear fit, resembling Fermi liquid behavior for the 
pressure above T c . The black solid line shows the pressure of 
a non-interacting Fermi gas, the dashed solid line the linear 
approximation, valid if T <C Tf- 

fit at high temperatures underestimates T c / \i c by 22%. 
This method resulted in T c /[i c — 0.32(2) in |S6| (shown 
for comparison), inconsistent with our determination. 
However, our value for T c / fx c — 0.40(3) agrees excel- 
lently with the most accurate theoretical determination 
T c /fi c = 0.400(14) from [SIS]. In a Fermi liquid, one has 

pfaT) = PoGM) (V 3/2 + ^7 1/2 ^ (^) 2 ) M- 

Although the agreement with the Landau Fermi liquid 
is fortuitous, we may nevertheless use it to model the 
data. From the linear fit versus (T/fj,) 2 , we obtain the 
fit parameters = 0.46(1), in agreement with the de- 
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termination of the same parameter from the chemical 
potential, and ^- = 1.04(2). Note that in the same tem- 
perature regime, a linear fit to the normalized pressure of 
a non-interacting Fermi gas (solid line in Fig. S4) gives 
^ = 0.91, while clearly, this value for a non-interacting 
Fermi gas should not differ from unity. The reason for 
discrepancy is that the pressure of the non-interacting 
Fermi gas in the temperature range 0.2 < (T//i) 2 < 1 
is not in the linear regime of low temperatures (shown 
as the dashed line in Fig. S4). Equivalently, the en- 
tropy per particle is not much smaller than ks- The fit- 
parameter m* should thus not be taken literally as the 
effective mass of particles in the normal state. However, 
the general finding of |S6) of an only weakly renormalized 
normal state at unitarity remains true. Compared to the 
non-interacting Fermi gas in the same regime of T/fj,, 
the unitary gas has an apparent effective mass enhance- 
ment of only about 12%. In view of the specific heat, 
which is not linear in temperature above T c , and the facts 
that neither T/n, nor the entropy per particle, are small 
above T c , we conclude that the usual Landau Fermi liq- 
uid picture, valid when the temperature is much smaller 
than the Fermi energy, cannot properly describe all ther- 
modynamic properties of the normal unitary Fermi gas. 
Note that our determination of fit parameters differs from 
the previous result [36] ^ = 1.13(3) and £„ = 0.51(2), 
possibly due to the mentioned differences in calibration 
between the experiments. This results in overall lower 
pressures (by a factor of about 13%) than in the present 
experiment, where we find £ = 0.376(5), to be shifted by 
at most 2% due to the error in the Feshbach resonance. 
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